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The magnetic quadrupole structure formation during the interaction of two ultra-short high power 
laser pulses with a collisionless plasma is demonstrated with 2.5-dimensional particle-in-cell simula¬ 
tions. The subsequent expansion of the quadrupole is accompanied by magnetic field annihilation 
in the ultrarelativistic regime, when the magnetic field can not be sustained by the plasma current. 

This results in a dominant contribution of the displacement current exciting a strong large scale 
electric field. This field leads to the conversion of magnetic energy into kinetic energy of accelerated 
electrons inside the thin current sheet. 
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Magnetic reconnection is a process which changes the 
topology of the magnetic field lines and allows the trans¬ 
fer of magnetic field energy into charged particles en¬ 
ergy m- It plays a fundamental role in the dynamics of 
magnetic confinement thermonuclear plasma [2] and has 
been considered as a plausible mechanism for high en¬ 
ergy charged particle generation in space plasmas [M]. 
This process is accompanied by the current sheet forma¬ 
tion [ 31 ], where the oppositely directed magnetic fields 
annihilate. Annihilation is a basic mechanism for recon¬ 
nection and has been investigated within the framework 
of dissipative magnetohydro dynamics [10]. In the limit of 
ultrarelativistic plasma dynamics the magnetic field an¬ 
nihilation acquires principally different properties related 
to the existence of a limiting value of the electron density 
m- This is due to the relativistic constraint on the up¬ 
per limit of the particle velocity which can not exceed the 
speed of light in vacuum and can sustain only a limiting 
magnetic field strength. 

The development of high power laser technology [T2| 
allows to access new regimes of magnetic field annihila¬ 
tion. When a high intensity laser pulse interacts with 
a plasma target the accelerated electron bunches gen¬ 
erate strong regular magnetic fields [ini HI]. The self¬ 
generated magnetic field in inhomogeneous near criti¬ 
cal density plasma also enhances fast ion generation[T5]. 
It has been predicted that nontrivial topology of self¬ 
generated magnetic field configurations should invari¬ 
ably lead to magnetic reconnection [T6] . Recent exper¬ 
iments have observed plasma outflows with MeV elec¬ 
trons and plasmoids generated in the reconnection cur¬ 
rent sheet [T71 [TS] . These experiments used long nano¬ 
second laser pulses and investigated reconnection in a 
parameter regime where the conditions for fast relativis¬ 
tic magnetic field line reconnection can not be reached. 

In this letter, a fast magnetic field annihilation in 


the relativistic collisionless regime driven by two paral¬ 
lel synchronized laser pulses with ultra-high intensity and 
femto-second pulse duration is studied with kinetic simu¬ 
lations. Different from Ref. m, we present the magnetic 
annihilation in the rear of the target where the plasma 
density is very low. The convection current is negligi¬ 
ble in the annihilation region and the variation of the 
magnetic field is compensated by the displacement cur¬ 
rent. The inductive electric field propagates forward and 
accelerates electrons and positrons in the current sheet. 

The Particle-In-Gell (PIG) simulations are performed 
with the relativistic electromagnetic code EPOGH m- 
Two s-polarized Gaussian pulses with peak intensity of 
10^^ W/cm^ are incident in the x direction and focused 
on the left edge of the target. By choosing s-polarization, 
the effects of the the high frequency laser field are mit¬ 
igated, which is helpful for observing clearly the self¬ 
generated magnetic field and the inductive electric field. 
The pulses durations are r = 15 fs and the laser spot 
sizes (EWHM) are 3 A, where A = 2'Kcju: = 1 lam is the 
laser wavelength and u is the laser frequency. The optical 
axes of the two laser pulses are at ^ = ±7 A. The sepa¬ 
ration guarantees the formation of two independent elec¬ 
tron bubbles, which do not overlap with each other. The 
simulation box is 195 Ax 340 A in x and y, respectively. 
The transverse size is large enough to avoid boundary 
effects. The preformed hydrogen plasma is located in 
20 A < X < 122 A with non-uniform density distribution 
in the x direction and uniform density in the y direction. 
The longitudinal density profile is shown in Eig. The 
plasma density linearly increases from 0 to 0.1 ric in the 
interval 20 A < x < 22 A, and then remains constant 
in 22 A < X < 62 A. Eor 62 X < x < 122 A, there is a 
downramp region, where the density linearly decreases to 
zero. Here ric = is the plasma critical density, 

rUe, e and Sq represent electron mass, electric charge, and 
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FIG. 1: (color online) Setup of the simulation. The longitudi¬ 
nal electric held values along y — ^ are taken at 33 To (green), 
69 To (red) and 144 To (blue), where To is the laser period. 
Zones I, II and III correspond to the regions of magnetic held 
generation, expansion and annihilation, respectively. 


vacuum permittivity, respectively. The transverse size of 
the target is 40 A. The mesh size is Sx = Sy = A/20. 
There are 10^ pseudo-particles in the simulation box and 
all the particles are initially cold. Open boundary condi¬ 
tions are employed for both particles and helds. 

The propagation of the laser pulses inside the plasma 
results in the formation of two electron bubbles due to the 
laser wakeheld effect ([21] and references cited therein). 
The strong wakehelds accelerate electrons in the longi¬ 
tudinal direction. According to Ampere’s law, the elec¬ 
tric currents produce magnetic helds, which is Bz in the 
2D case. As a result, a magnetic quadrupole conhgura- 
tion is formed as shown in Fig. [^a). At this moment, 
the two magnetic dipoles do not touch each other in the 
vicinity of the central axis (^ = 0). The maximum mag¬ 
netic held can be calculated by using Faraday’s law as 
V X B = /ioJe + fioSod'E/dt with /io being the vacuum 
permeability. Assuming the quasistatic conditions, we 
take the magnetic dipole transverse size to be of the or¬ 
der of R ^ Cy/^jujp^ which corresponds to the radius of 
the self-focusing hlament for the laser pulse with the nor¬ 
malized amplitude ao = eE/rrieUJC^ where the Langmuir 
frequency is uOp = {Anne^ . Then the magnetic 

held can be estimated as Bz = iiQU^ec^/ujp with the 
Lorentz factor 7 « ao [E]. It gives 5;^ ^ 1.8 x 10^ T in 
our case, which is consistent with the simulation result. 

When the bubbles propagate into the density down- 
ramp region, where n{x) = —no{x — 122 A)/(60 A), both 
the bubble size and the magnetic held dipoles expand 
transversely since they experience forces acting on the 
vortex proportional to Vn x Fi [^. Here ft is the poten¬ 
tial vorticity. Figure 2(b) shows the prohle of Bz along 
X = 27 A, 40 A and 68 A at 33 To, 51 To and 75 To, 
respectively. During the time between 33 To and 51 To, 
the quadrupole is still propagating in the uniform den¬ 
sity region and the displacement between the blue solid 
peak and the green dashed peak is about 0.3 A. Dur- 



FIG. 2: (color online) (a) The z-component of the magnetic 
held (Bz) distribution at t = 63 Tq. (b) The expansion of 
the magnetic dipole structures. The blue solid line, green 
dashed line and red dotted line represent the Bz prohle along 
X = 27 A, 40 A and 68 A at 33 To, 51 To and 75 To, respectively. 


ing this stage the growth of the magnetic held amplitude 
(about up to 8500 T) is signihcant since more and more 
electrons are trapped and accelerated by the wakehelds. 
The situation changes, when the quadrupole enters the 
inhomogeneous region, where the displacement between 
the green dashed peak and the red dotted peak becomes 
about 1.6 A. At the same time, the magnetic held maxi¬ 
mum amplitude slightly decreases which is in accordance 
with Ertel’s theorem [23]. The shift of the peak posi¬ 
tions proves that the magnetic helds are expanding in 
the transverse direction. As a result, the two magnetic 
dipoles approach each other near the central axis, creat¬ 
ing steep gradients in the magnetic helds which in turn 
facilitate annihilation. 

The central axis becomes the line where the opposite 
magnetic held lines annihilate and rearrange the topol¬ 
ogy. In the vicinity of the central line, the current sheet 
is formed, which is along the x direction. Magnetic anni¬ 
hilation in the current sheet region is accompanied by the 
generation of the inductive electric held. It is regarded 
as an important signature of magnetic reconnection. The 
longitudinal electric held distribution along x = 65 A, 
80 A and 90 A lines at the instants of t = 81 Tq, 99 Tq and 
123 To respectively are plotted in Fig. [^a). At 81 Tq and 
99 To, the intense electric held is located near the laser 
axes {y = ±7 A). Here the electric held corresponds to 
the plasma wave: Eq = cmeOOpje 10 GV/cm. However, 
the longitudinal electric held in the current sheet (near 
(^ = 0)) becomes much higher at t = 123 To, while the 
magnitude of the electric held along laser axes still keeps 
the same order. The extra contribution comes from the 
magnetic annihilation which releases the magnetic energy 
converting it into electric helds. The magnetic held dis¬ 
tribution at 138 To is shown in Fig. 3(b). We note that 
a region in the current sheet around 100A<x<110A, 
where the magnetic quadrupole breaks, has unique prop¬ 
erties. Behind the breaking region, the magnetic held is 
quite smooth with continuous distribution. However, in 
front of the breaking region, the magnetic held becomes 
hlamented and disrupted, which indicates that the mag¬ 
netic held lines with opposite directions are annihilating. 
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FIG. 3: (color online) (a) The growth of the inductive longi¬ 
tudinal electric field. The blue solid line, green dashed line 
and red dotted line represent the Ex profile along x = 65 A, 
80 A and 90 A at t = 81 To, 99 To and 123 To, respec¬ 
tively. (b) Bz distribution at 138 To. (c) The profiles of 
Ex inside the current sheet at 132 To (blue), 138 To (green) 
and 144 To (red), respectively, (d) Contributions of dif¬ 
ferent terms in the Faraday’s law at 138 To along the x- 
direction and the spatial averaged profile inside the current 
sheet (—A < y < X)\ -^ < {V x 'B) >x (blue), the convection 
current density < Je >cc (green) and the displacement current 
density < Jd >cc (red). Inset: The smoothed value near the 
region where the inductive electric field grows. 


As a result, this filamented region also corresponds to the 
region of location of strong inductive longitudinal electric 
field. This is well seen in the E^ profiles in Fig. 3(c). The 
green line shows the same moment with (b) and has the 
peak {Ex > 40 GV/cm) around 110 A < x < 115 A. Fig. 
3(c) also shows the evolution of the longitudinal electric 
field in the current sheet. The strong inductive electric 
field moves forward almost with the speed of light. 

To investigate the regime of the inductive electric field 
growth, we compare the contributions of different terms 
in Faraday’s law. The profiles of (I/to) < (V x B) >a^, 
the convection electric current density < Je >cc= —enVx 
and the displacement current < Jd >a:= s^dtEx are 
shown in Fig. 3(d). In the region 80 A < x < 110 A, 
there are return electrons which induce a strong convec¬ 
tion current. Here (V x B)^^ is balanced by < Je >cc. 
However, in the region oi x > 110 A, the electron density 
is very low due to the downramp distribution. The vari¬ 
ation of magnetic fields can no longer be compensated by 
the convection current. Therefore, (V x B)^, is balanced 
by the displacement current. This induces the growth of 
the inductive electric field. The high frequency oscilla¬ 
tions of < Jd >x also agree well with the filaments in 
magnetic field distribution. After smoothing the oscilla¬ 
tion (see Fig. 3(d), inset), (V X B)^ becomes precisely 
in accordance with the displacement current within the 
region of the strong inductive electric field. 
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FIG. 4: (color online) The distributions of the electron longi¬ 
tudinal momentum px (averaged over one cell) at 69 To and 
138 To are plotted in (a) and (b), respectively. 


The electrons within the current sheet get extra energy 
compared to the electrons outside of the current layer due 
to the contribution of the inductive electric field. Figures 
4(a) and (b) show the distributions of electron longitudi¬ 
nal momentum px at 69 Tq and 138 Tq, which represent 
the situation without and with magnetic annihilation ef¬ 
fects, respectively. The momenta of the electrons within 
the current sheet {—2 X < y < 2 X) and that of the elec¬ 
trons near the outer wings of laser axis (8 A < ^ < 25 A 
and —25 X < y < —8 A) are comparable in magnitude at 
t = 69 To. The contributions to Px at this moment come 
from the electron oscillations in the plasma waves and 
from the return electrons along the bubble shell. How¬ 
ever, at t = 138 To with the process of magnetic annihi¬ 
lation, a strong backward accelerated electron bunch is 
formed in the current sheet. The maximum momentum 
reaches —110 rUgC, while the momentum growth of the 
electrons in the wings is not significant. Furthermore, 
the origin of the backward accelerated electrons in the 
current corresponding to the action of the strong dis¬ 
placement current, i.e. to the inductive electric field as 
is shown in Figure 3. This also demonstrates that the 
electron acceleration is driven by the magnetic annihila¬ 
tion. Such a kind of violent electron acceleration in the 
current sheet region is a clear evidence for the magnetic 
annihilation. 

The inductive electric field can also be used to accel¬ 
erate efficiently positrons. We perform a simulation with 
a mixture plasma of 10% of positrons and 90% of pro¬ 
tons. Positron generation by laser plasma interaction has 
been demonstrated in Ref. [24] . The separation between 
the two laser pulses is enlarged from 14 A to 20 A. If 
the pulse separation is too small, all the positrons en¬ 
countered by the laser pulses are pushed forward by the 
ponderomotive force. As a result, when the inductive 
electric field grows, there are no positrons left in the cur¬ 
rent sheet. When the pulse separation is enlarged, a frac¬ 
tion of the positrons in the current sheet is not pushed 
forward. They are squeezed together by the transverse 
ponderomotive force and form a high density region be¬ 
tween the two bubbles. With the growth of the inductive 
electric field, they are accelerated forward up to high en¬ 
ergy and forms a high density positron bullet. All the 
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FIG. 5: (color online) The density and momentum x — px 
distribution of positrons at 171 To are plotted in (a) and (b), 
respectively. Here Upo = 0.01 ric is the initial positron density, 
(c) The energy spectrum of the inductive held accelerated 
positron bullet at 120 To (blue), 150 To (green) and 180 To 
(red). 


other simulation parameters are the same as in the sim¬ 
ulations discussed above. Figure 5(a) shows the positron 
density distribution at 171 Tq. The high density peak 
around x = 170 X indicates positrons directly acceler¬ 
ated by the laser ponderomotive force. Behind that, the 
high energy positron bullet is accelerated by the inductive 
electric held. The bullet has maximum momentum equal 
to Px ~ 350 rUeC, which is much higher than the positron 
momentum in the hrst peak corresponding to the laser 
direct accelerated positrons (see Fig. 5(b)). The energy 
spectrum evolution of the positron bullet is plotted in 
Fig. 5(c). The peak energy increases from 20 MeV to 
160 MeV in 60 Tq. As the spectrum shows, the positron 
bullet is quasi-monoenergetic with the energy spread of 
about SEjE « 13%. The angular divergence approxi¬ 
mately equals 14 °. The bullet contains the charge of 
44 pC//im. The energy spectrum shows that the in¬ 
ductive field accelerates the positrons to higher energy 
than is the case for the electrons. This is because the 
inductive electric field moves forward with the propagat¬ 
ing laser pulses continuously accelerating positrons. By 
contrast, the accelerated electrons move backwards and 
quickly leave the region of the inductive field. 

In conclusion, this letter identifies a new regime of colli¬ 
sionless relativistic magnetic annihilation using petawatt 
lasers operating on very fast time scales. In the magnetic 
annihilation region, the variation of V x B can only be 
compensated by the displacement current. The associ¬ 
ated electric field in the current sheet accelerates elec¬ 
trons to very high energy, opposite to the laser propaga¬ 
tion direction. Similarly, the mechanism can efficiently 
accelerate positrons in the forward direction, which could 
be used as a signature of this new regime in an experi¬ 
ment. This paper provides predictive simulations for the 
upcoming petawatt laser installation like ELI [25] . 
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